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\q , Abstract. In this paper we calculate the Galactic positron flux from dark matter 

annihilation in the frame of supersymmetry, taking the enhancement of the flux by 
existence of dark matter substructures into account. The propagation of positrons in 
the Galactic magnetic field is solved in a realistic numerical model GALPROP. The 
£SJ ' secondary positron flux is recalculated in the GLAPROP model. The total positron 

■ flux from secondary products and dark matter annihilation can fit the HEAT data well 

when taking a cuspy density profile of the substrctures. 
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+3 . 1. Introduction 



The astronomical observations indicate that most of the matter in our universe is dark 
(see for example PQ). The evidences come mainly from the gravitational effects of 



the dark matter component, such as the rotation curves of spiral galaxies [2J, the 
gravitational lensing [3] and the dynamics of galaxy clusters [4]. The studies of 
primordial nucleosynthesis [5] , structure formation [6J and cosmic microwave background 
(CMB) [7] show that this so-called dark matter (DM) is mostly non-baryonic. The 
nature of the non-baryonic dark matter remains one of the most outstanding puzzles in 
particle physics and cosmology. 

A large amount of theoretical models of DM have been proposed in literatures [8J. 
All these candidates of non-baryonic DM particles require physics beyond the standard 
model of particle physics. Among the large amount of candidates, the most attractive 
scenario involves the weakly interacting massive particles (WIMPs). In particular, the 
minimal supersymmetric extension of the standard model (MSSM) provides an excellent 
WIMP candidate as the lightest supersymmetric particle (LSP), usually the neutralino, 
which are stable due to R-parity conservation [I]. 

WIMPs are possible to be detected beyond the gravitational effects. One viable 
way is to detect the elastic scattering signals of DM particles on the detector nuclei 
(direct searches) [9j [101 [TTJ [12]. This is the most sensitive method at present for 
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DM detection |13j . Another way is to search the self-annihilation prodcuts of the DM 
particles (indirect searches), such as neutrinos [H], 7-rays [151 US], antiprotons [17] and 
positrons [HI [19]. Since the positrons, antiprotons and diffuse 7-rays are secondary 
particles in the cosmic rays with low fluxes, these kinds of DM annihilation products 
are easier to be distinguished from the astrophysical background. Direct and indirect 
detections of DM particles are complementary ways to each other. 

The balloon-borne instrument High-Energy Antimatter Telescope (HEAT) was 
designed to determine the positron fraction over a wide energy range with high statitical 
and systematic accuracy. Based on the three flights(1994, 1995 and 2000), this 
collaboration reported an excess of positrons with energies higher than several GeV and 
peaked in the range 7 ~ 10 GeV [20J [21]. It seems that the astrophysical production 
cannot give enough positrons as observed [22] . and there might exist exotic sources 
of positrons. It has been pointed out that the excess may indicate the signal of DM 
annihilation [T8l [19] . However, Baltz et al showed that for supersymmetry (SUSY) DM 
a "boost factor" > 30 has to be introduced to fit the HEAT data [23]. It was thought 
that the local dumpiness of DM distribution might account for this "boost factor". 
However, Hooper et al argued that it was less possible for a DM clump to be close 
enough to the Earth to contribute sufficient positron flux [23]. 

Stimulated by the recent numerical simulation which shows that the lightest 
subhalos in the Galaxy could be extended to 10~ 6 M Q with a huge amount of about 10 15 
[25] , we try to calculate the Galactic positron flux again in the frame of supersymmetry 
taking the enhancement by the large amount of subhalos into account. There was effort 
try to explain the HEAT data by the nearby mini-clumps [26]. However, the possibility 
of this scenario is also extremely low [27J. Another difference of this work from the 
previous studies most of which adopt the Green's function method to calculate the 
positron's propagation is that we calculate its propagation in the Galactic magnetic 
field using the numerical package GALPROP [28J , which adopts realistic distribution of 
interstellar gas and radiation field. Our result shows that the two methods lead to some 
difference. To give a consistent result we also calculate the secondary positron flux using 
GALPROP. The propagation parameters are adjusted to fit observations such as B/C, 
proton and electron spectra. Our results show that the secondary positron has better 
agreement with the HEAT data than which adopted in previous studies [T8], [19j [23], [26] 
and decreases the tension between predictions and data. 

The paper is orgnized as following: in Sec. [2] we present the SUSY model. In Sec. 
[3] we give our treatment of subhalos distribution and density profiles. Sec. HI introduces 
the positron's propagation and the GALPROP model. The results are given in Sec. [H 

2. Positron Production from Neutralino Annihilation 

The LSP, generally the lightest neutralino, is stable in the R-parity conservative MSSM 
and provides a natural candidate of dark matter. It is a combination state of the 
super symmetric partners of the photon, Z° boson and neutral Higgs bosons. 
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Positrons can be produced in several neutralino annihilation modes. Most of them 
come from the decay of gauge bosons produced in chanels XX ~^ ZZ, XX ~~ > W + W~ , or 
from the cascades of final particles, such as fermions and Higgs bosons. The spectrum 
of positrons depends on the neutralino mass and its annihilation modes. There is also a 
direct channel XX e + e~ and produces monochrome positrons with energy E e + = m x . 
However, the branching ratio of this channel is usually small. In the following discussion 
we neglect this "line" contribution to the positron spectrum. 

The source function of positrons from DM annihilation can be written as 



where a is the positron generating cross-section, dn/dE is the positrons spectrum in 
one annihilation by a pair of neutralinos and p(r) is neutralino density distribution in 
space. The source term Q(E e +,r) is given in unit of GeV~ 1 m~ 3 s~~ 1 . 

The source term is calculated in MSSM by doing a random scan using the software 
package DarkSUSY In these SUSY models we choose a few models which satisfy 
all the experimental bounds and give large fluxes and appropriate spectrum so that we 
can fit the HEAT data well. However, there are more than one hundred free SUSY 
breaking parameters even for the R-parity conservative MSSM. A general practice is to 
assume some relations between the parameters and greatly reduce the number of free 
parameters. For the processes related with dark matter production and annihilation, 
only seven parameters are relevant under some simplifying assumptions, i.e., the higgsino 
mass parameter fi, the wino mass parameter M 2 , the mass of the CP-odd Higgs boson 
rriA, the ratio of the Higgs Vacuum expectation values tan/3, the scalar fermion mass 
parameter mj, the trilinear soft breaking parameter A t and Af,. To determine the low 
energy spectrum of the SUSY particles and coupling vertices, the following assumptions 
have been made: all the sfermions have common soft-breaking mass parameters mr, all 
trilinear parameters are zero except those of the third family; the bino, wino, and gluino 
have the mass relations, M\ = 5/3 tan 2 9 W M 2 , M 3 = a 3 (M^)/a e m sin 2 9\yM 2 , coming 
from the unification of the gaugino mass at the grand unification scale. 

The parameters are constrained in the following ranges: 50GeV < 

M 2 , M a , nif < lOTeV; 1.1 < tan/3 < 55; -3m 5 - < A t , A b < 3niq; sign(p) = ±1. 
The SUSY models are required to satisfy the theoretical consistency requirement, such 
as the correct vacuum breaking pattern, the neutralino being the LSP and so on. The 
accelerator data constrains the models further from the spectrum requirement, the in- 
visible Z-boson width, the branching ratio of b — > 57 and the muon magnetic moment. 

The constraint from cosmology is also taken into account by requiring the relic 
abundance of neutralino 0.091 < Q x h 2 < 0.118, which corresponds to the 3a bound from 
the cosmological observations [7j. The effect of coannihilation between the fermions is 
taken into account when calculating the relic density numerically. 




(1) 
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3. Dark Matter Distribution 



To determine the positron source term from DM annihilation in Eq. (OQ) we have to 
specify the DM density profile p(r). Based on the N-body simulation results the DM 
density profile can be written in a general form as [30] 

P= (r/r s )^[l + (r/r^O 3 ^)/"' ^ 
where p s and r s are the scale density and scale radius respectively. These two parameters 
are determined by the measurement of the virial mass of the halo and the concentration 
parameter determined by simulation. The NFW profile was first proposed by Navarro, 
Frenk and White (31] with (a, (3, 7) = (1, 3, 1). However, Moore et al. gave another form 
of DM profile with (a, (3,j) = (1-5,3, 1.5) to fit their simulation result [32]. The Moore 
profile has steeper slope near the Galactic center than the NFW profile. There are also 
simulations showing that the density profile may not be universial [331 El]- Reed et 
al. show that 7 = 1.4 — 0.08 log(M/M*) increases for halos with smaller masses [34J. 
Therefore we also adopt a 7 = 1.7 profile for the whole range of DM halos for simplicity. 
At present there are controversials on which one represents the actual density profile. 



100 



z = 
median 

-1 sigma 
+1 sigma 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

1e-05 1 100000 1e+10 1e+15 

M vir( h ~ lM SUN> 

Figure 1. c v j r as a function of subhalo mass m sv fo predicted by the Bullock model 
35j. The median and the la values are given. 

For a DM halo with mass M v j r , we adopt a semi-analytic model of Bullock et al [35j . 
which describes the concentration parameter function of viral mass and redshift, 

to determine the parameters in the density profile. We plot the Cyj r — M v j r relation in 
Figure [H with the la range. The median c v j r — M v j r relation at redshift zero is adopted. 
The scale radius is then determined as 7g w = r v ir/ c vir' r™ OOTe = r§ w /0.63 and r§ = 



The Galactic positron flux and dark matter substructures 



„nfw 



/(2 — 7). The virial radius r v j r in above formula is defined by M v j r = 47rAp c r^ r /3 
with A = 200 and the critical density of the universe pc = 139M //cpc 3 . Once r s is 
fixed, the scale density p s can be got by solving the equation M v j r = / p(r)d 3 r. In 
Figure [2] we show the scale radius r s , scale density p s , virial radius r v j r , surface mass 

density p(r v j r ), effective volume £ = J(p/ p ) 2 d 3 r, intrinsic boost factor B c = j-^^^ 
as functions of masses M v ; r , for the DM halos with NFW mass distributions. We adopt 



a reference mass density po — 0.3Ge V/cm 3 . Each of these charicteristic quantities is 
scaled by its maximum value in the mass range 2 x 10 -6 ~ 2 x 10 10 M Q , the minimal and 
maximal masses of the Milky Way (MW) subhalos we adopt in the following discussion. 
We list these scale factors as: r max = 2.72kpc, pg 11 ^ = 3.52 x 10 9 M Q /kpc 3 , r^ x = 
55.62fcpc, p(r vir ) max = 1.52 x 10~ 4 GeV/cm 3 , £ max = 1909.11A;pc 3 , 5 max = 39.89. 
The behaviors of these quantities can be understood as following. For NFW profile, 



an analytic calculation shows that r v j r oc and r s oc M^/cy^.. Because the 

concentration parameter c v j r is weakly dependent on the mass according to Figured], rs 
is approximately proportional to M^?. The scale density ps is nearly oc c 3 ^./ ln(l+c v j r ), 



rl/3 



which shows a slight decrease as the increase of halo mass, while the surface density 
P( r vir) ~ P( s )l ' c vir ^ Vl n (l + c vir) ^ s a l mos t constant. The effective volume £ is 



approximately scaled with M, 
little smaller for higher mass halos, because it is less concentrated for higher mass halo. 



l vir c vir' 



and the intrinsic boost factor B c oc c^. B c is a 
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Figure 2. Some characteristic quantities of DM halos with NFW mass distribution 
as functions of masses. 



Having known the general properties of the DM halos given above, we now 
turn to the discussion of the MW subhalos. High resolution simulations have 
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revealed that a large number of self-bound substructures survived in the galactic halos 
[361 EZl EH], [391 SQl SH S21 S3] . Especially a recent simulation conducted by Diemand et 
al. [25] showed that the first generation objects as light as the Earth mass can survive 
until today. The number of such minihalos is huge, reaching 10 15 . The existence of a 
wealth of subhalos will greatly enhance the DM annihilation flux, since the annihilation 
rate is propotional to the density square as shown in Eq. ([1]) . 

The simulations show that the number density of subhalos can be approximately 
given by an isothermal profile with a core [S] 

n(r) = 2n H (l + (r/r H ) 2 )- 1 , (3) 

where njj is the relative number density at the scale radius r#, with r# being about 0.14 
times the halo virial radius r# = 0.14r v j r . The result given above agrees well with that 
in another recent simulation by Gao et al. [45] . Eq. ([3]) shows that the radial distribution 
of substructures is generally shallower than the density profile of the smooth background 
in Eq. ([2]) . This is due to the tidal disruption of substructures which is most effective 
near the galactic center. 

The mass function is well fitted to the simulation result as dn/dm gub oc m ^^- We 
then get the number density of substructures with mass ^ sub at the position r to the 
galactic center 

n(m sub , r) = n (^-) (1 + {r/r H fy l , (4) 

where M v j r is the virial mass of the MW, n is the normalization factor determined by 
requiring the number of subhalos with mass larger than 1O 8 M is about 500 in a halo 
with M = 2 x 10 12 M Q [38]. 

Taking the subhalos into account the density square in Eq. (CQ) is given by 

p 2 (r)^(p 2 (r))=^ mooth (r) + (p 2 ub (r)), (5) 

where (Pg ub ( r )) means the average density square of subhalos according to the 
distribution probability. Since there is no correlation between the mass and spatial 
distribution in Eq. (J4j) we get (Pg ub ( r )) a ^ the position r is given by an integral of mass 

/" m max / r \ 

^sub( r )) = J ' n Kub> r ) (J Plnb dV ) - dm sub 
m min 

flQ /•'"max 

= 1 + (r/r H y L n ^Kub)"d-sub > (6) 

where refers to density of the subhalo at r and V is its volume. We defined a 

function /(^ sub ) in the equation above as 

^subr 1 - 9 r Kub) .^, 



/KubMlf^ L 47rr'^ ub (r')dr', (7) 



M vir J Jo 

which depends only on the subhalo mass since the virial radius and the density profile are 



all determined by the subhalo mass. The minimal subhalos can be as light as 10 6 M 



[25] while the maximal mass of substructures is taken to be 0.01M v j r 
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We notice that there are unphysical singularities at the center of the subhalo 
which may lead to /( m SUD ) divergent for the index 7 equal to or greater than 1.5. 
A cutoff r cu ^ is introduced within which the DM density is kept a constant value 
Pmax due to the balance between the annihilation rate and the rate to fill the 
region by infalling DM particles [17] . Applying typical setting of parameters, we get 
Pmax ~ 10 18 — 10 19 M G /kpc 3 . In the following discussion, pmax = 10 18 M Q /kpc 3 is 
adopted for NFW and Moore profiles while we consider pmax = 10 18 — 2 x 10 19 M Q /kpc 3 
for the profile with 7 = 1.7. 

From Eq. ® we know that, (p 2 ub ) oc / /(m sub ) dm gub = / ™ su b / ( m S ub ) d lnm sub- 
The quantity rn sn ^ x /(m su |,) is plotted in Figure EJ which shows the relative contribu- 
tion to the annihilation signals from subhalos at different mass bins at the logarithmic 
scale. We can see from this figure that subhalos of different mass bins contribute almost 
the same to the annihilation products. By extending the minimal subhalo mass from 
the present numerical resolution of about 10 6 M Q to 10~ 6 M Q , the contribution will be 
~ 4 times larger. 
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Figure 3. x /(^ig^) as a function of subhalo mass rrig^. 

Figure H] shows (p 2 ) 1 ^ 2 for different DM profiles. The existence of subhalos will 
increase (p 2 ) 1//2 significantly at large radii. Because the annihilation source is density 
square dependent, the further dumpiness of DM will effectively increase the annihilation 
signals, and can play the role of a "boost factor". The more concentrated of the halo 
profile, the larger of the (p 2 ) 1//2 value. While at small radii there is no enhancement 
compared to the smooth component for the tidal disruption makes the subhalos be 
difficult to survive in inner Galaxy. We can see that the model with the most cuspy 
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Figure 4. (p 2 (r)) 1 / 2 distribution of Galaxy for various DM profiles. 

subhalos adopted here has ~ 3 times larger (p 2 ) 1 ^ 2 than the smooth one at the solar 
location, which means that the positron flux observed on the Earth will be higher by 
about an order of magnitude. 

4. The Positron Propagation 



A complexity in calculating the positron flux is that positrons are scattered by the 
Galactic magnetic field (GMF) and we have to calculate its propagation in GMF. 
The propagation of charged particles in the GMF is diffusive. In addition they will 
also experience energy loss processes due to ionization and Coulomb interaction in 
the interstellar medium; for electrons and positrons there are additional synchrotron 
radiation in the GMF, bremsstrahlung radiation in the interstellar medium and inverse 
Compton scattering in the interstellar radiation field. The interaction between cosmic- 
ray particles and the interstellar medium (mostly Hydrogen and Helium) will lead to 
the fragmentation of nuclei. For radioactive nuclei the decay should be also taken into 
account. 

The full propagation equation including the convection of cosmic-ray particles and 
reacceleration processes is written in the form 



d_ 

dp 



pip 



op op p z 

Tf T r ' 



(8) 



where ip is the density of cosmic ray particles per unit momentum interval, Q(x.,p) is 
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the source term, D xx is the spatial diffusion coefficient, V is the convection velocity. 
The third term describes the reacceleration using diffusion in momentum space, p is the 
momentum loss rate, 77 and r r are timescales for fragmentation and radioactive decay 
respectively. 

The interaction between cosmic-ray particles and the interstellar medium produces 
the secondary particles, such as positrons, antiprotons and some other heavy elements. 
The source term of the secondary cosmic-rays is given by 

Q(x,p) = (3cipp(x,p)[a H (p)n H (^) + o- He {p)n H e{^)\, (9) 

where t/> p (x, p) is the density of primary cosmic-rays, (3c is the velocity of particles, o~h 
and o~H e are the production cross sections for the secondary particles from the progenitor 
on H and He targets, nn and nu e are the interstellar Hydrogen and Helium number 
densities, respectively. 

The propagation equation (JHI) can be solved analytically under some simplification 
assumptions [HJ US]. A numerical solution to Eq. (JSJ is given in the GALPROP 
model developed by Strong and Moskalenko [28] , which takes all the relevant processes 
into account. The realistic distributions for the interstellar gas and radiation fields are 
adopted in GALPROP. The detail of this model can be found in |28j . 

The secondary-primary ratio such as B/C depends on the propagation process. 
Therefore the propagation parameters are adjusted to describe the B/C ratio, the 
electron and proton spectra et al . We include the nuclei up to Z = 28 and relevant 
isotopes. The injection spectra of protons and heavier nuclei are assumed to have 
the same power-law form in rigidity. The nuclei injection index below and above the 
reference rigidity at 9 GV are taken as 1.98 and 2.42 respectively. The electron injection 
index below and above 4 GV are taken as 1.60 and 2.54 respectively. For propagation, 
we use the diffusion reacceleration model [49J. The spatial diffusion coefficient is taken 
as pD (p/p ) 5 , where D = 5.4 x 10 28 cm s" 1 , p = 4 GV, and 5 = 0.33. The Alfven 
speed is va = 30 km s -1 . The height of the propagation halo is taken as Zh = 4 kpc. 

In Figure. [5] we show the expected B/C from the propagation model. The 
propagation model describes the observational data very well. We will use this model 
to calculate the secondary positron spectrum from cosmic rays fragmentation given in 
Eq. (Q and the primary positron propagation from DM annihilation given in Eq. ([1]). 

5. Results and Discussion 

Once the SUSY model is specified and the contribution from subhalos is taken into 
account we can incorporate the source term of primary positrons in GALPROP and 
calculate its spectrum on the Earth. In Figure [6] we show the propagated positron 
fluxes for three different SUSY models with m x = 121, 195, 242 GeV. For comparison 
the local source terms E 2 Q(E, r = r ) according to Eq. ([I]) are also shown. We see that 
propagation makes the spectra some different from the source spectra. This is because 
the propagation is energy dependent. The source spectra are somewhat different among 
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Figure 5. Calculated Boron-to-Carbon ratio compared with observational data. The 
lower line represents the local interstellar value, while the upper one is the result after 
solar modulation with $ = 400MV. Observational data: ACE [SO] . Ulysses [51] . 
Voyager 52J, HEAO-3 [53], for others please refer to [54 . 



the three models because these models lead to different annihilation final states. In the 
calculation a Moore profile for the smooth component and 7 = 1.7 profile for subhalos 
with p max = 2 x 10 19 M Q /kpc 3 are adopted. 

In Figure [7] we show the positron flux for the m x = YllGeV model for different 
DM profiles. It indicates that the fluxes can be different for several times among various 
DM profiles. The cuspier profiles lead to higher positron fluxes. It should be noted that 
because of the propagation effects the spectra from different DM profiles are different, 
although they have the same source spectrum. At lower energy the fluxes have bigger 
differences, since the low energy positrons have longer propagation distance, and may 
trace the (p 2 ) farther away from the Earth. The differences between the NFW (or 
Moore) profile with and without subhalos are not significant. For subhalos with 7 = 1.7 
profile the flux is greatly enhanced. 

We then compared the results between GALPROP and other propagation methods. 
The results are shown in Figure El The fluxes for m x = 121 and 242 GeV based 
on the Green's functions given by Kamionkowski & Turner [55] and Baltz & Edsjo 
[T8] are plotted. The BE Green's function is to solve the diffusive propagation 
equation analytically and similar to the GALPROP method. However, GALPROP uses 
realistic astrophysical inputs, and can contain more complex physical process such as 
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Figure 6. The propagated positron fluxes on the Earth for three neutralino models. 
Moore profile, containing subhalos with 7 = 1.7 and p max = 2 x f0 19 M©/fcpc 3 is 
adopted. Upper lines show the local annihilation sources E 2 Q(E, r) in unit of GeV 
to _3 s _1 for the three neutralino models. 
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Figure 7. Positron flux on Earth for various DM profiles. The mass of neutralino is 
121GeV. 
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Figure 8. Comparison of positron fluxes of different propagation models. 

reacceleration, while in BE's method only an average energy loss rate is considered and 
without convection or reacceleration. After adjusting the propagation parameters to be 
equivalent, we get similar results between these two models. KT's result shows a larger 
difference from these two models. This is comprehensible that the KT Green's function 
is acquired by solving the leaky-box equation of propagation with a position-independent 
source. At low energies the discrepancy is larger because of longer propagation path. 
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Figure 9. Positron fraction e + / (e + + e~) calculated in this work. The left panel shows 
the result of a Moore DM profile without substructures. The right one is the same as 
the left but including the subhalos with 7 = 1.7 and p max = 2 x 10 19 -M©/fcpc 3 . Data 
points are HEAT measurements in the three fights in 1994, 1995 and 2000 [20| [21]. 



To give a consistent Galactic positron flux we calculate the background (secondary) 
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positron fraction generated by interaction between primary cosmic-rays and interstellar 
gas in the same propagation model using GALPROP. In the Figures [9] and [10] we give 
the background positron fraction. Our results give better description to HEAT data 
compared with the adapted background fraction in [181 EE 1231 [26] . Therefore the new 
results requries a smaller "boost factor" to the DM annihilation signals. 

Adding the positron flux coming from DM annihilation to the secondary positrons 
we get the total positron flux. The positron fraction e + /(e + + e~) is plotted in Fig. [9] 
for m x = 121 GeV. We can see from this figure that the expectation reproduces the 
measurements relatively well if DM subhalos are taken into account. The model with 
subhalos can contribute several times more than which without subhalos. The best fit 
to the HEAT data requries "adjustment factor" of 9.2 and 2.5 for signals from DM 
annihilation without and with subhalos respectively. 

We don't think the discrepancy of a factor 2.5 is serious since there are large 
uncertainties in calculating the positron fluxes. Firstly there are uncertainties for 
the positron propagation, as shown in Figure El Secondly the random distribution 
of DM subhalos may have large variance as discussed in [27] since positrons lose energy 
quickly and have large fluctuations. The variance may lead to large positron flux due 
to accidence, while in this work we just give the average positron flux. Thirdly we can 
choose SUSY model which produces more positrons, as shown in Figure [10] for a model 
with m x = 195GeV. The best fit to HEAT data requries a factor of 1.2 for this model. 
It is remarkable that without introducing a large "boost factor" nor the nonthermal 
production of neutralinos as done in [19] we could explain the HEAT data quite well. 
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Figure 10. The same as Fig. [9] but for neutralino with mass m x = 195GeV. 
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It is interesting to notice that at about 8GeV, HEAT measurements show some 
fine structures of the positron fraction, which cannot appear in our calculation. Maybe 
there are some specific sources or some specific physical process to generate positrons 
at this energy. However, the errors of the data are too big to give a definite assertion of 
this property. If it is confirmed by further experiments such as AMS-2, there should be 
deeper studies about the models we used here. 

In summary, in this work we give a consitent and detailed calculation of the Galactic 
positron flux. We recalculate the background positron flux from cosmic-ray collisions 
with the interstellar gas in a realistic propagation model GALPROP. The results soft 
the discrepancy between data and expectations. For the primary positrons from DM 
annihilation we take the enhancement of subhalos into account. The result shows that 
the HEAT data can be explained by the two components without introducing "boost 
factor" or nonthermal production of DM. To account for the excess of positron flux 
extending subhalo masses to 1O~ 6 M and requiring a cuspy profile as 7 = 1.7 are 
necessary. 
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